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ABSTRACT 


The  Phase  1  study  of  the  "Ravleiph"  critical  angle  phenomenon 
has  been  completed.  An  automated,  computer  controlled  goniometer 
was  constructed  and  perfected.  Controlled  measurements  on  a  variety 
of  materials,  including  materials  of  interest  to  the  NOT  community, 
were  carried  out  and  compared  with  theoretical  predictions  made  by 
a  number  of  authors. 

Experimental  results  indicate  that  these  theories  are  only 
first  approximations  to  reality.  We  find  a  number  of  significant 
experimental  anomalies,  such  as  the  appearance  of  more  than  one 
amplitude  dip.  Certain  theoretical  questions  can  also  he  raised 
regarding  these  theories  including  the  question  of  momentum 
conservation  at  the  interface.  Both  the  experimental  and  theoretical 
observations  suggest  that  we  are  not  yet  in  possession  of  a  complete 
theory . 

Phase  II  plans  are  designed  to  investigate  both  experimental lv 
and  theoretically  these  questions  so  that  our  goal  of  a  better 
theoretical  understanding  of  this  phenomenon  can  be  achieved.  Given 
this  improved  understanding,  practical  quantitative  NOT  techniques 
based  on  the  "Rayleigh"  critical  angle  phenomenon  will  become 


available. 
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INTRODUCTION  AND  BACKGROUND 


The  introduction  of  high  performance  machines,  such  as  aircraft 
engines  and  gas  turbines,  demands  increasingly  exotic  materials 
capable  of  handling  high  temperature  and  stress.  These  materials  are 
generally  of  the  ceramic  or  super  alloy  type  formed  by  powder  metal¬ 
lurgy.  Since  these  materials  are  intended  to  operate  under  severe 
environments,  it  is  important  to  be  able  to  characterize  them  as  much 
as  possible.  Such  parameters  as  grain  size,  residual  stress,  and 
number  and  size  of  defects  have  an  important  bearing  on  the  machin- 
abilitv  and  performance  of  any  component  fabricated  from  the  material. 
Thus,  there  is  a  need  for  methods  to  perform  a  nondestructive 
evaluation  of  these  materials. 

Current  methods  for  nondestructive  evaluation  of  materials 
include  x-ray,  eddy  current,  dye  penetrant,  and  ultrasonics.  All  of 
these  techniques  have  advantages  and  disadvantages.  One  of  the  most 
popular  methods  is  ultrasonic,  used  in  the  pulse-echo  mode.  This 
method  sends  pulses  of  ultrasound  into  the  material  and  examines  the 
echos  returned  from  within.  The  size  of  internal  defects  is  estimated 
by  analyzing  the  echos.  It  is  also  possible  to  create  images  by 
scanning  the  transducer  over  the  material  surface,  and  recording  the 
echos  in  their  correct  spatial  position. 

Ultrasonic  inspection  has  several  disadvantages  when  applied 
to  the  new  materials.  One  of  the  major  problems  is  the  inability  to 
resolve  near-surface  defects  due  to  ringing  and  reverberation 
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phenomena*:.  The  other  major  problem  is  a  lack  of  resolution  (on 
the  order  of  1  mm).  For  critical  components,  near-surface  defects 
are  the  most  important  since  fractures  occurring  at  the  surface  will 
cause  rapid  deterioration  and  can  lead  to  catastrophic  failure  at 
the  high  temperatures  and  pressures  which  must  be  endured.  In 
addition,  current  ultrasonic  methods  are  incapable  of  measuring  grain 
size  and  surface  residual  stress  which  also  influence  material 
behavior . 

The  purpose  of  the  research  reported  here  is  to  investigate 
the  feasibility  of  harnessing  a  nearly  unused  phenomenon  of  ultrasound, 
the  so-called  Rayleigh  critical  angle  effect,  to  provide  measurement 
of  grain  size,  defect  location  and  size,  and  surface  residual  stress 
distribution  to  name  a  few.  If  feasibility  can  be  established,  the 
way  is  open  to  the  development  of  inspection  instrumentation  capable 
of  evaluation  materials  as  they  come  from  the  production  line. 

It  is  well  known  that  an  ultrasonic  wave  propagating  from  one 
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medium  to  another  undergoes  significant  alteration  at  the  boundary. 

When  a  wave  propagates  from  a  liquid  into  a  solid,  three  new  waves 
generally  appear;  one  reflected  and  two  propagated.  The  propagated 
waves  travel  at  angles  defined  by  Snell's  law,  and  the  reflected  wave 

u 

obeys  the  law  of  reflection.  Specifically,  the  direction  of  the 
propagated  bulk  longitudinal  wave  is  given  by: 

sin  0 1  ■>  =  sin  On  ,  (1) 

vi  i 


-2- 


the  propagated  bulk  shear  wave  angle  by: 

sin  022  =  sin  Oji  ,  (2) 

vi  1 

and  the  reflected  wave  angle  by: 

0ir  =  0  1 1  O) 

where  vn  =  velocity  of  longitudinal  waves  in  liquid, 

v  1 2  =  velocity  of  longitudinal  waves  in  solid, 

v  2  2  ~  velocity  of  shear  waves  in  solid,  and 

Oil  =  angle  of  incidence. 

The  above  equations  hold  for  ideal  media  having  no  attenuation. 
When  attenuation  is  present,  Snell's  law  becomes  more  complicated 
since  velocity  and/or  wave  number  now  becomes  a  complex  quantity. 

Snell's  law  predicts  the  interesting  phenomenon  known  as 
the  critical  angle.  This  is  easily  seen  from  an  examination  of 
equations  1  and  2.  At  some  incident  angle,  the  angle  of  refraction 
of  either  the  longitudinal  or  shear  wave  becomes  90°  indicating 
that  one  or  both  of  these  refracted  waves  no  longer  exist  in  the 
material.  Typically,  for  steel,  the  longitudinal  critical  angle 
occurs  at  about  17°,  and  the  shear  critical  angle  occurs  at  about 
31°.  The  simple  equations  predict  no  more  than  this.  In  fact, 
the  classical  theory  for  the  behavior  of  sound  predicts  complete 
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reflection  at  and  bevond  the  shear  critical  angle. 

Experiments  show,  however,  that  the  behavior  just  beyond 
the  shear  critical  angle  is  much  different  than  predicted  by 

f-  8 

classical  theory.  Ihe  amplitude  of  the  reflected  wave  dips 
sharply  just  bevond  the  shear  critical  angle  and  its  phase 
undergoes  a  relatively  large  shift,  typically  360°.  It  is  observed 
that  just  beyond  (approximately  1°)  the  shear  critical  angle  a 
Rayleigh-tvpe  surface  wave  is  set  up  on  the  liquid-solid  boundary. 
Obviously,  the  theory  must  be  modified  to  account  for  this  behavior. 

3-  i  6 

This  has  been  done  by  several  authors  who  have  introduced 

attenuation  in  this  problem.  This  results  in  a  new  reflection 
coefficient. 

The  models  seem  to  predict  observed  behavior  for  unbounded 
beams  with  good  accuracy.11  Both  the  dip  in  amplitude  and  the 
phase  change  seem  to  be  predicted  quantitatively  as  well  as 
qualitatively  in  cases  where  it  has  been  experimentally  tested. 

Such  unbounded  beam  theories  have  also  been  extended  to  include  the 
case  of  bounded  beams  which  leads  to  further  insights  into 
experimental  observations. 

In  addition  to  describing  the  behavior  of  the  reflected  wave 
at  the  critical  angle,  these  theories  indicate  that  there  should 
exist  a  frequency  at  which  there  is  zero  reflection  at  the  critical 
angle,  called  the  frequency  of  least  reflection  (FIR).  This  point  is 

apparently  dependent  upon  grain  size,  and  appears  to  present  a 

J  3 

method  for  its  measurement.  The  theory  of  Richardson  also  seems 


to  predict  that  grain  size  could  be  measured  to  a  resolution  on 

1  f. 

the  order  of  0.2  pm  by  measuring  the  FLR  to  0.1%. 

In  the  course  of  the  experimental  verification,  some  interesting 
transient  effects  were  observed  by  several  authors.1  It  was  noted 
that  the  phase  response  near  the  critical  angle  involved  a  time  lag. 
That  is,  it  took  a  finite  amount  of  time  after  initiation  of  the 
ultrasonic  burst  before  the  phase  of  the  reflected  wave  changed  from 
the  classically  predicted  value  to  that  obtained  with  the  steady 
state  theories.  It  is  speculated  that  this  time  delay  may  be  a 
measure  of  the  viscious  properties  of  the  material.  Thus,  it  appears 
that  extension  of  these  theories  to  include  transient  effects  might 
be  profitable. 

Any  changes  that  affect  attenuation  and  velocity  can  be 
detected  by  the  behavior  of  the  reflected  wave  at  the  critical  angle. 
Some  general  statements  can  be  made  about  these  effects. 

1)  The  depth  of  the  reflectivity  minimum  depends  on 
attenuation,  principally  that  of  the  shear  wave  although 
the  phase  shift  is  insensitive  to  changes  in  shear  wave 
attenuation. 

2)  Angular  position  of  the  surface  wave  critical  angle 
depends  on  velocity,  principally  that  of  the  shear 
wave . 

3)  The  two  effects  listed  above  appear  to  be  independent. 

4)  Changes  in  critical  angle  with  frequency  probably  result 
from  velocity  gradients  with  depth. 

The  critical  angle  phenomenon  is,  however,  a  near  surface 
effect  extending  into  the  material  to  a  distance  of  roughly  the 
wavelength  of  surface  waves.  Thus,  the  selection  of  frequency 


could  be  used  to  probe  different  near-surf ice  layers.  For 
example,  it  may  be  possible  to  measure  the  stress  distribution  in 
the  solid  as  a  function  of  depth  by  tracking  the  critical  angle  as 
a  function  of  frequency.  Martin  has  shown  that  the  surface  wave 
velocity  can  be  measured  by  monitoring  the  critical  angle  with 
sufficient  accuracy  to  enable  discrimination  of  3.3  x  105  dyn/cm7 
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uniaxial  stress.  By  starting  at  a  high  frequency  and  progressing 
to  lower  frequencies,  the  stress  variation  with  depth  could  be 
analyzed  layer  by  layer. 

There  are  a  variety  of  ways  in  which  the  critical  angle 
phenomenon  could  be  used  to  produce  images.  Among  these  are: 

A.  Plot  of  critical  angle  versus  poistion, 

B.  Plot  of  the  FLR  at  the  critical  angle  versus  position, 

C.  Plot  of  amplitude  of  the  reflected  wave  for  fixed  angle 
of  incidence  at  the  nominal  critical  angle  versus 
position,  and 

D.  Plot  of  phase  of  the  reflected  wave  for  fixed  angle  of 
incidence  at  the  nominal  critical  angle  versus  position. 

The  various  imaging  modes  listed  above  (A,  B,  C,  D)  would  be 
sensitive  to  different  combinations  of  material  properties  as 
indicated  below: 

Mode  A:  Material  "variation''  and/or  surface  stress 

Mode  B:  Grain  size 

Mode  C:  Grain  size,  surface  stress,  material  variation,  and 

Mode  D:  Material  variation  and/or  surface  stress 
Material  variation  may  consist  of  changes  in  hardness,  presence  of 


defects,  changes  in  chemical  composition,  etc.;  in  general, 
anything  that  causes  local  variation  of  the  shear  wave  velocity. 

Only  one  known  experiment  in  critical  angle  imaging  has  been 
2  0 

reported.  In  this  experiment,  method  D  was  used  to  image  a 
hardness  variation  and  a  subsurface  defect  with  excellent  results. 
The  possibilities  discussed  above  are  only  the  obvious  ones  based 
on  the  results  available  at  this  time.  Other  applications  will 
arise  as  more  work  is  done. 

In  the  present  report,  the  results  of  both  experimental  and 
theoretical  studies  of  this  problem  are  presented.  Certain 
anomalous  results  and  various  theoretical  difficulties  suggest  that 
the  correct  theoretical  description  is  not  yet  available.  Studies 
are  proposed  to  remedy  these  deficiencies. 
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2.0 


REVIEW  OF  PHASE  I  PI.ANS 


There  were  five  principal  tasks  outlined  in  the  original 
proposal  for  the  first  vear's  work. 


2 . 1  Year  I  Tasks  As  Planned 

In  brief  these  tasks  were: 


Task  1  -  Theoretical  Review  of  Steady  State  Theory 

X  3  i  i, 

The  theories  of  Richardson,  Merkulova,  and  Bertoni  and 
Tamir,?1-'3  will  he  reviewed  and  combined  to  provide  the  most 
complete  picture  of  the  critical  angle  phenomenon  possible. 

As  well,  a  search  of  the  modern  literature  will  be  made  to 
ascertain  that  no  new  theories  have  been  proposed  since 
their  work.  The  Russian  literature,  in  particular,  will 
be  carefully  perused  since  much  theoretical  work  is  going 
on  there.  At  the  end  of  this  task,  the  steadv  state 
theory  should  he  firmly  established  and  ready  for  exploita¬ 
tion. 


Task  2  -  Experimental  Backup  for  Task  1 

A  precision  goniometer  will  be  designed  and  fabricated  so 
that  precise  critical  angle  measurements  may  be  made.  A 
system  consisting  of  the  goniometer,  transducers,  ultrasonic 
generator,  and  measurement  instruments  such  as  a  vector 
voltmeter  will  be  assembled.  With  the  advice  of  the  Materials 
Research  Center  of  the  NBS,  a  number  of  flat  samples  of  well 
characterized  materials,  such  as  Ab  6061-T6,  will  be  obtained. 
These  will  be  completely  characterized  (with  respect  to  the 
critical  angle  parameters)  by  recording  the  phase  and  amplitude 
of  the  reflected  wave  versus  angle  with  frequency  as  a 
parameter.  From  these  curves  other  curves  mav  be  derived. 

These  are  critical  angle  versus  frequency,  amplitude  minimum 
at  the  critical  angle  versus  frequency,  width  of  amplitude 
dip  around  the  critical  angle  versus  frequency,  and  slope 
of  phase  change  at  the  critical  angle  versus  frequency. 

From  these  derived  curves,  certain  characteristics  of  the 
material  can  be  established.  These  are  amplitude  and 
frequency  of  least  reflection,  frequency  of  narrowest  ampli¬ 
tude  dip,  and  frequency  of  largest  phase  slope.  It  is  expected 
that  the  frequencies  should  all  be  the  same. 
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Tho  ultrasonic  parameters  of  the  sample,  velocity,  and 
attenuation  of  the  longitudinal  and  shear  waves  would 
be  inserted  into  the  theory  and  the  same  theoretical 
curves  drawn. 

Thus,  Task  2  would  provide  a  thorough  comparison  of  theory 
and  practice. 

T ask_  3  -  Development  of  Tr ansient  Theory 

The  theory  reviewed  and  revised  in  Task  1  would  be  extended 
to  include  the  tr  insient  case.  It  is  suspected  that  the 
time  lag  to  steady  state  behavior  mav  be  an  additional  and 
independent  variable  useful  in  characterizing  the  material. 

Task  4  -  Experimental  Backup  for  Task  3 

The  same  samples  and  measurement  system  developed  in  Task  2 
would  be  used  to  verify  the  transient  theory.  Instead  of 
continuous  wave  signals,  step  functions  (gated  sine  waves) 
would  be  used  with  careful  attention  paid  to  the  time  interval 
before  steady  state  response  is  achieved.  Since  at  this  time 
it  is  not  known  exactly  what  to  expect,  a  particular  set  of 
experiments  cannot  be  specified. 

Task  5  -  Applications  Experiments 

Since  it  is  already  established  that  subsurface  flaws  can 
be  detected  by  the  phase  shift  at  the  critical  angle,  this 
application  will  be  investigated  first.  A  sample,  similar 
to  those  used  in  Tasks  2  and  4,  will  he  prepared  by  drilling 
flat-bottomed  holes  !nto  it.  Those  holes  will  be  arranged 
in  specific  patterns  of  size,  distance  from  the  surface,  and 
sepurat  i on .  ihe  goniometer  will  be  mounted  on  an  x-y  scanner 
so  that  a  two-dimensional  distribution  of  critical  angle 
parameters  mav  be  detected  and  displayed.  In  this  wav, 
resolution  and  sensitivity  of  tho  various  parameters  to  the 
t law  can  be  ascert a  ined . 
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3.0 


PROGRESS  TOWARD  ACHIEVING  GOALS 


In  the  course  of  making  acoustic  critical  angle  measurements 
early  in  the  program,  it  gradually  became  obvious  that  a  better 
experimental  apparatus  was  needed.  Great  difficulty  was  experienced 
in  leveling  and  orienting  samples.  Sample  surfaces  were  often 
found  to  be  too  rough  or  insufficiently  flat  to  obtain  reproducible 
measurements  and  it  was  also  determined  that  results  are  temperature 
dependent.  Measurements  made  at  different  times  were  often  not 
repeatable  owing  to  unknown  temperatures. 

These  difficulties  have  been  overcome  by  incorporating  a 
gimbaled  sample  holder  (see  Figure  1),  and  employing  a  temperature 
controlled  bath.  We  have  also  obtained  a  grinding  lap  for  producing 
flat  samples  with  uniform  surface  characteristics. 

3 . 1  Description  of  Apparatus 

A  gimbaled  sample  holder  (Figure  1)  was  designed  and  constructed. 
This  device,  which  is  patterned  after  a  microscope  mirror  mount,  is 
designed  to  accept  samples  measuring  not  more  than  two  inches  in 
maximum  dimension  and  possessing  at  least  one  flat  surface. 

The  sample  holder  is  fitted  with  a  threaded  ring  and  various 
disc  "windows"  for  adapting  to  samples  of  different  sizes  and 
shapes.  Not  only  can  the  orientation  of  the  plane  sample  surface  be 
changed,  but  the  elevation  may  also  be  changed  by  slight  adjustments 
of  the  rings.  With  this  new  device  in  place,  most  of  the  former 
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FIGURE  1 


THE  EXPERIMENTAL  SYSTEM  FOR  MAKING  MEASUREMENTS  OF  R(0) 
AND  4)  (0 )  VS  P  CONSISTS  OF  A  GONIOMETER  IN  FIGURE  1(a), 
AND  A  G  1MBALF.D  SAMPLE  HOLDER  IN  FIGURE  1(b).  THE 
COMPUTER  CONTROLLED  AUTOMATED  SYSTEM  IS  SHOWN  IN  FIGURE 
1(c)  AND  1(d).  DATA  IS  STORED  ON  DISC  FOR  ARCHIVAL 
PURPOSES  AND  MAY  HE  DISPLAYED  AS  SHOWN  IN  FIGURE  1(c). 
THE  SYSTEM  IS  COMPATIBLE  WITH  PDP  11/23,  APPLE,  OR  IBM 
PERSONAL  COMPUTERS.  WHEN  x-v  SCANNING  IS  PERFORMED, 

THE  G 1MBALED  SAMPLE  HOLDER  IS  REPLACED  BY  A  SAMPLE  MOUNT 
THAT  CAN  BE  SCANNED.  SCHEMATIC  DIAGRAMS  SHOWING  THE 
APPARATUS  AND  SUPPORT  ELECTRONICS  ARE  PROVIDED  IN  FIGURE 
2.  THE  SYSTEM  IS  IMMERSED  IN  A  SMALL  WATER  TANK  (NOT 
SHOWN)  TO  ACHIEVE  ACOUSTIC  COUPLING  TO  THE  SAMPLE. 


( 


1(a)  GONIOMETER  AND  SAMPLE  HOLDER 
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problems  with  sample  leveling  and  orient  ini;  have  apparent  1  v  been 
el im inn  ted  . 

In  addition  to  the  gimbaled  sample  holder,  several  other 
minor  mechanical  improvements  were  made.  The  goniometer  apparatus 
was  f  itted  with  lei;  mounts  (see  Figure  1),  and  a  mtu  h  smaller  water 
tank  was  constructed.  This  reduced  water  volume  allows  the  system 
to  he  placed  on  a  lab  bench,  and  also  reduces  problems  that  would 
be  associated  with  maintaining  a  large  volume  of  water  at  a  fixed 
t  emperatur e . 

1.2  Computer  Hardware 

One  of  tiie  principal  requirements  of  Phase  1  is  to  use  a 
device  of  the  foregoing  type  to  obtain  reliable  data  for  assessing 
the  validity  of  various  theories  of  the  critical  angle  phenomenon. 
Measurements  of  reflected  amplitude  R(h)  and  phase  .*(  ')  versus  the 
angle  of  incidence  0  are  required. 

We  have  been  using  the  PDF  11/2)  computer  system  to  make  and 
store  measured  data  using  our  automated  measuring  system.  This 
computer  is  unnecessarily  large  for  this  task.  Accordingly,  we 
have  made  plans  for  adapting  the  system  to  a  smaller  computer. 

In  Phase  II  we  plan  to  employ  an  IBM  personal  computer  instead  of 
the  POP  11/23. 

) .  3  Set_  Up  Procedures 

With  the  new  mechanical  improvements  in  place,  we  have  been 


able  to  develop  a  simple  set  up  procedure  that  allows  one  to  make 
rapid  and  reliable  measurements.  Several  tests  of  repr oduc ib i 1 i t v 
involving  the  same  and  different  observers  and  the  same  and  different 
samples  have  been  made.  There  are  still  effects,  including  possible 
nonlinear  effects  to  be  discussed,  that  mav  at  times  affect  reli¬ 
ability.  However,  we  feel  that  we  now  have  the  ability  to  make 
measurements  that  are  sufficiently  reliable  so  that  we  can  begin  to 
make  meaningful  comparisons  between  theorv  and  experiment. 

3.4  The  Automated  Measuring  System 

The  original  proposal  called  for  the  construction  (Task  2) 
of  a  goniometer  of  unspecified  design  for  making  acoustic  critical 
angle  measurements.  Not  only  was  such  a  device  constructed  and 
the  necessary  auxiliary  support  equipment  assembled  (see  Figures  1 
and  2),  but  it  has  been  converted  to  a  computer-controlled  automated 
measuring  system.  Ive  may  program  the  system  to  begin  measurements 
at  some  initial  angle  of  incidence  0^  and  continue  in  steps  of  AG 
until  a  final  angle  9.  is  reached. 

Data  is  automatically  analyzed  to  find  the  angle  corresponding 
to  the  minimum  reflection  amplitude  (the  critical  angle),  and  is 
displayed  graphically  (see  Figure  3)  on  a  CRT  along  with  descriptive 
data  regarding  the  experiment  and  the  sample.  The  data  is  also 
stored  on  disc  for  archival  purposes.  The  system  allows  us  to  make 
very  accurate  measurements  of  R(0)  and  <J)(0)  in  a  rapid  fashion. 

The  computer-controlled  system  is  also  capable  of  operating  in  an 
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FIGURE  2  (Continued).  SCHEMATIC  DIAGRAM  SHOWING  APPARATUS  USED  TO  MEASURE  THE 

R EFFECTED  AMPLITUDE  R(<->  AND  THE  PHASE  I  (")  OF  THE  REFLECTED 
WAVE.  THE  SYSTEM  ALLOWS  ONE  TO  PERFORM  AUTOMATED  COMPUTER 
CONTROLLED  MEASUREMENTS  AND  DISPLAY  THEM  ON  A  CRT  MONITOR. 


x-y  scanning  mode:  thereby,  making  possible  "images"  of  the  scattered 
wave  at  some  fixed  angle,  usually  the  critical  angle. 

In  Figure  4  we  illustrate  an  example  of  the  capability  of 
our  device  to  produce  x-y  scanned  images  of  R(6)  and/or  $(0)  at  some 
fixed  angle  0  (usually  a  critical  angle).  This  capability  (a  Task  5 
goal)  has  also  been  used  on  holes  at  varying  depths;  and  it  has  been 
shown  that  by  varying  the  frequency,  the  depth  of  penetration  of  the 
surface  wave  can  be  varied  to  determine  the  depths  to  such  holes. 

Part  of  Task  2  involved  attempting  to  determine  material 
parameters  by  fitting  experimental  data  to  the  Richardson  model. 
Accordingly,  we  have  acquired  a  complete  set  of  programs  (courtesy 
of  R.  L.  Richardson)  which  allows  one  to  calculate  the  reflection 
coefficient  R(0)  and  phase  <J>(0)  versus  0  where  0  is  the  angle  of 
incidence  of  an  unbounded  longitudinal  wave  on  a  liquid-solid 
interface. 

These  programs  were  rewritten  for  use  on  our  PDP  11/23 
computer  system  and  curves,  such  as  that  shown  in  Figure  5  for  a 
stainless  steel  water  interface,  were  computed  and  plotted. 

Various  input  parameters  such  as;  material  attenuation,  shear  wave 
velocity,  and  so  on  Wire  chosen  for  these  plots. 

We  hope  eventually  to  be  able  to  fit  such  curves  to  our 
experimental  data;  and  thereby,  infer  the  material  parameters. 
Standing  in  the  way  of  this  are  two  problems.  First,  as  mentioned 
elsewhere,  we  are  not  yet  apparently  in  possession  of  a  complete 
theory.  Second,  our  method  of  measuring  data  involves  a  complex 
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FIGURE  3.  TYPICAL  EXPERIMENTAL  CURVES  FOR  ALUMINA.  R(0) 

IS  THE  REFLECTED  AMPLITUDE  VS  INCIDENT  ANGLE  Q. 
$(0)  IS  THE  PHASE  OF  THE  REFLECTED  WAVE.  THE 
CRITICAL  ANCLE  FOR  SURFACE  WAVE  GENERATION  OCCURS 
HERE  AT  ABOUT  32°. 
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'  FIGURE  4.  A  FLAT  BOTTOMED  1/4  INCH  DIAMETER  HOLE  WITHIN  1/32 

INCH  OF  THE  SURFACE  IN  ALUMINUM  IS  IMAGED  IN  THE 
x-y  SCANNING  MODE  (SEE  FIGURE  2).  AT  THE  CRITICAL 
4  ANGLE  =  Pt.r,  THE  HOLE  IS  VISIBLE  IN  BOTH  PHASE 

-1(0)  AND  AMPLITUDE  R(P)  PLOTS.  THIS  ILLUSTRATES 
THAT  A  SURFACE  WAVE  IS  GENERATED  AT  0  =  "cr . 
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bounded  beam  (spherical  focused  beam)  rather  than  an  unbounded  beam 
(Richardson's  programs  cover  this  case). 

Nevertheless,  for  many  materials,  curves  of  R(0)  and  <f>(0) 
versus  G  look  qualitatively  similar  to  those  of  Figure  5.  Furthermore, 
although  we  use  a  complex  bounded  beam,  we  also  employ  a  very 
narrow  angle  "point"  detector  which  should  see  only  one  reflected 
wavevector  (unless  many  wavevectors  interact  nonlinear ly) .  All  beam 
shift  effects  apparently  occur  at  points  sufficiently  distant  from 
our  detector;  that  is,  this  single  detected  wavevector  may  well  act  like 
a  reflected  unbounded  beam  and  should  be  comparable  in  many  cases 
with  Richardson's  calculations.  We  will  accordingly  attempt  to 
curve  fit  several  well  behaved  experimental  materials  to  the 
Richardson  model;  and  thereby,  fix  the  material  parameters  used  in 
his  model. 

3. 5  Literature  Search 

The  goal  of  the  Task  1  literature  search  was  to  study  and 
compare  works  on  bounded  and  unbounded  acoustic  beams  as  they 
influence  the  observed  critical  angle  phenomenon. 

Accordingly,  the  University  of  Washington  information 
retrieval  service  among  other  sources  was  used  to  obtain  references 
pertaining  to  the  Rayleigh  critical  angle  phenomenon  and  its 
applications.  A  search  of  three  data  bases  was  carried  out.  These 
were:  A)  COMPENDEX  -  January  1970-present,  935,688  records.  Sources 
include  3500  journals,  publications  of  engineering  societies  and 
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organizations,  papers  from  proceedings  of  conferences  and  selected 
government  reports,  and  books;  B)  INSPEC  -  1969-present,  1,661,133 
citations.  Sources  include  Physics  Abstracts,  Electrical  and 
Electronic  Abstracts,  and  Computer  and  Control  Abstracts;  C)  NTIS  - 
1964-present,  833,513  citations.  Sources  include  government 
sponsored  research,  development,  engineering  plus  analyses  prepared 
by  federal  agencies,  their  contractors  or  grantees.  Agencies  such 
as  NASA,  DD,  DOE,  HEW,  HUD,  AND  DOT  are  included.  In  addition,  we 
have  obtained  a  number  of  useful  references  from  other  published 
papers  and  from  conversations  with  investigators  in  the  same  field. 

The  existing  literature  on  this  phenomenon  tends  to  fall 
into  two  categories:  A)  Theories  of  unbounded  incident  plane  waves, 
concerned  primarily  with  the  development  of  material  models  which 
explain  the  observed  amplitude  dips  and  phase  shifts;  and  B)  Theories, 
where  the  incident  waves  are  bounded,  which  attempt  to  describe  what 
happens  when  a  bounded  rather  than  an  unbounded  beam  is  involved.  In 
principle,  all  theories  of  this  type  are  extensions  of  the  unbounded 
theories.  Their  primary  goal  is  to  explain  a  variety  of  puzzling 
nonspecular  effects,  such  as  the  beam  shift  effect  first  observed 
by  Schoch . 6 

In  section  5  we  will  review  two  of  the  more  interesting  and 
successful  theories,  namely  that  of  Richardson  and  Becker13  and  that 
of  Pitts^4’  and  similar  theories. 26-2 9 
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4.0 


SUMMARY  OF  RESEARCH  RESULTS  (PHASE  I) 


In  the  following  section,  we  describe  experimental  results 
or.  a  set  of  controlled  and  uncontrolled  samples.  General  statements 
regarding  the  qualitative  features  of  the  data  are  made  before 
discussing  the  prospects  of  direct  comparison  of  these  results  with 
existing  theory. 

4 . 1  Samples 

Existing  literature  on  the  subject  of  the  critical  angle 
scattering  phenomenon  has  focused  largely  on  a  theoretical  explanation 
with  only  modest  experimental  support.  That  is  not  to  say  that 
published  experiments  do  not  agree  with  published  theories;  rather 
the  number  and  types  of  materials  examined  does  not  seem  to 
constitute  a  very  serious  test  of  these  theories.  Accordingly,  we 
undertook  a  preliminarv  survey  of  the  experimental  properties  of 
a  wider  class  of  materials.  We  have  examined: 

1)  Naturally  occurring  minerals  and  rocks  (uncontrolled) 
including  single  crystals  of  quartz,  fluorite,  and 
calcite 

2)  Single  crystals  (controlled)  of  silicon  and  saphire 

3)  Alloys  (controlled  and  uncontrolled)  of  aluminum,  iron, 
titanium,  copper,  and  brass 

4)  Composites  such  as  an  aluminum  carbon  fiber  composite 

5)  Glasses  such  as  fused  silica  and  obsidian 

6)  "Ceramics"  such  as  alumina,  tungsten  carbide,  and 
tungsten  powder  composites 
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Only  the  alloys,  ceramics,  and  composites  are  of  direct 


interest  to  the  NDT  community;  however,  our  understanding  of  the 
theory  can  be  improved  by  examining  the  response  of  simple  single 
crystals  rather  than  complex  polycrystals  such  as  alloys.  We 
included  rocks  and  minerals  in  this  survey  in  order  to  discover 
possible  anomalous  behavior  due  to  the  complex  structure  and  chemical 
composition  often  exhibited  by  these  naturally  occurring  materials. 

4. 2  Temperature  Control 

All  experiments  were  performed  at  a  temperature  of  25°C. 

This  was  essential  for  accurate,  reproducible  measurements  as 
illustrated  in  Figure  6  where  the  critical  angle  for  a  titanium  sample 
is  plotted  vs.  temperature  for  T  =  5  to  25°  C. 

Because  the  apparatus  may  exhibit  some  distortions  with  a 
change  of  temperature  and  because  the  surface  and  bulk  velocities 
(of  the  sample  and  the  water  bath)  are  functions  of  temperature, 
there  are  considerable  changes  in  0  with  changes  in  T.  In  Figure  6 
for  Ti,  we  find  A0/AT  =  (.1°)/(°C).  The  heat  bath  employed  was 
sufficiently  stable  to  allow  measurements  of  0  to  about  *.025°  or 
better  as  shown  in  Figure  7. 

4.3  Exp e rime ntal  Procedure 

Samples  with  existing  flat  but  rough  surfaces  were  prepared  by 
grinding  these  surfaces  with  # 100  and  #600  carborundum  in  that  order. 
The  sample  was  then  placed  in  the  gimbaled  sample  holder  and  the 
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Critical  Angle  (Deg) 


FIGURE  6 


1 


0c  vs  T(°0  at  3.1  MHz 
IT i  Sample) 


TEMPERATURE  EFFECTS.  A  SAMPLE  OF  TITANIUM  WAS  EXAMINED 
AT  A  SERIES  OF  TEMPERATURES  RANGING  FROM  5  TO  25°C. 
STRAINS  IN  THE  APPARATUS  AND  CHANGES  IN  VELOCITY  IN 
BOTH  THE  WATER  BATH  AND  THE  SPECIMEN  CONTRIBUTE  TO 
CHANGES  IN  THE  CRITICAL  ANGLE  0C.  IT  IS  CLEAR  THAT 
TEMPERATURE  VARIATIONS  ARE  IMPORTANT  IN  SUCH  MEASUREMENTS 
SINCE  WE  OBSERVE  A  CHANGE  OF  APPROXIMATELY  1°  IN  0C 
FOR  5°  C  (BETWEEN  10-15°  C)  CHANGE  IN  TEMPERATURE. 
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Critical  Angle  (Deg) 


FIGURE  7 


V 


0c  vs  TIflE  (7075  AL,  3.8  rife,  75  D 

0C 


0  1  2 


Time  (Hours) 


REPRODUCIBILITY  OF  MEASUREMENTS.  AN  ALUMINUM  SAMPLE 
WAS  MAINTAINED  AT  25°  C  AND  MEASUREMENTS  OF  0C  WERE 
MADE  OVER  A  PERIOD  OF  HOURS.  A  NUMBER  OF  MEASUREMENTS 
WERE  ALSO  MADE  AFTER  REPEATEDLY  REPLACING  AND  REMOVING 
THE  SAMPLE.  THESE  MEASUREMENTS  (NOT  SHOWN)  HAD  ESSENTIALLY 
THE  SAME  ERROR  AS  THE  SAMPLE  THAT  WAS  NOT  REMOVED  DURING 
THE  SERIES  OF  MEASUREMENTS. 
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,;on  iometer  placed  in  a  .‘mall  water  tank.  The  water  (clean  tap  water) 
is  circulate..!  in  a  temperature  control  reservoir,  and  all  measurements 
are  made  at  When  fresh  water  is  initially  placed  in  the  reservoir, 

circulation  for  several  hours  is  required  to  eliminate  dissolved  gasses 
that  would  result  in  bubble  formation  on  the  transducers  and  the 
sample  surface. 

When  the  sample  is  in  place  and  the  system  temperature  has 
equilibrated,  the  goniometer  is  calibrated  by  setting  up  on  a  known 
angle  ( l''  =  10°)  which  is  the  initial  state  assumed  bv  the  computer- 
control  software. 

After  the  appropriate  electronic  settings  are  made,  the 
.•inhaled  sample  holder  is  properly  oriented.  This  is  done  by  instructing 
the  goniometer  from  the  computer  keyboard  to  move  to  an  incident  angle 
b  =  45°.  This  angle  is  usually  far  from  any  critical  angle  for  most 
materials  (see  Table  I).  If  it  does  happen  to  fall  on  a  critical 
angle,  one  can  set  up  at  any  other  angle  provided  6  4  O^wheic 

is  the  surface  wave  critical  angle, 
e  r 

The  gimbaled  sample  holder  is  then  roughly  leveled  by  eye. 

Next,  fine  adjustments  are  made  on  the  gimbaled  sample  holder  until 
tiie  reflected  signal  amplitude  R(0)  is  a  local  maximum  at  the  chosen 
angle.  Wtien  this  occurs,  the  sample  is  assumed  to  be  properlv 
oriented  in  space  and  a  measurement  can  then  be  made. 

A  reconnaissance  run  on  the  sample  is  then  made  to  locate 
the  critical  angle.  Data  is  plotted  on  a  CRT  and  examined  visually 
to  locate  the  angle.  More  appropriate  initial,  final  and  angle 
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increments  can  then  be  keyed  in  if  necessary  and  a  final  run  made. 

The  entire  procedure,  given  a  properly  prepared  sample,  takes 
less  than  five  minutes  per  sample.  Data  stored  on  disc  cartridges 
can  be  called  up  for  analysis  or  examination  at  a  later  date. 


4 . 4  Qualitative  Features  of  Experimental  Data 

Before  discussing  what  may  be  called  "anomalous"  experimental 
results,  let  us  describe  the  general  qualitative  features  of  those 
sample  measurements  that  are  more  or  less  similar  to  theoretical 
predictions.  In  Figure  8,  we  illustrate  schematically  a  "typical" 
curve  for  R(0)  and  4>(9)  vs.  the  incident  angle  0. 

At  point  A,  0  is  equal  to  the  shear  critical  angle;  that  is, 
the  angle  0S  at  which  the  shear  wave  propagates  along  the  water-solid 
boundarv.  From  Snell's  law,  this  condition  reads 


sin  0S 
v 


sin  tt/2 


(4) 


w 


0  =  sin  1  (v  /v  ) 

s  w  s 


(5) 


where  v  is  the  longitudinal  velocity  in  water  and  v  is  the  shear 

W  ’  s 


wave  velocity  in  the  solid. 

From  A  to  B ,  the  amplitude  R(0)  begins  to  drop  and  a  shift  in 
the  phase  <}>(0)  begins.  At  point  B,  R(0)  approaches  a  minimum  and  the 
rate  of  change  of  the  phase  is  maximal.  Finally,  at  point  C  there  is 
a  local  maximum  in  R(ft)  before  R(0)  settles  down  to  a  more  or  less 
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These  results  are  similar  but  not  identical  to  theoretical 

13,  15.  2U 

predictions  for  unbounded  beams,  It  is  well  to  keep  in 

mind  the  fact  that  we  are  examining  a  very  small  portion  of  a 
nonuniform  reflected  beam  that  results  from  a  special  kind  of 
bounded  beam  (a  focused  beam).  Figure  9  illustrates  the  experi¬ 
mental  arrangement  we  employed  which  involves  a  focused  lens  system 
and  a  "point"  receiver  (.1  mm).  This  point  receiver  is  designed  to 
examine  only  one  "wavevector"  out  of  the  range  of  wavevectors  incident 
on  the  sample.  As  such,  this  experimental  arrangement  to  a  certain 
extent  eliminates  the  problem  of  examining  the  beam  profile.  Only 
the  component  that  is  in  the  location  of  the  "would  be"  specular 
reflection  along  the  axis  of  the  lens  is  being  examined.  Trailing 
waves,  which  result  from  the  reradiation  of  the  surface  wave  back  into 
the  water,  are  largely  eliminated. 

Basically,  the  point  receiver  should  produce  results  comparable 
to  an  unbounded  beam.  Therefore,  comparisons  between  our  results  and 
the  unbounded  beam  theories  appear  to  be  appropriate,  at  least  in  a 
first  approximation. 

4 . 5  Results  and  Anomalies 

A  rather  "typical"  experimental  curve  is  illustrated  in 
Figure  10.  The  specimen  was  a  titanium  alloy  (6AL-2Cb-lTa-lMO) ,  and 
measurements  were  performed  at  25°C  and  3.1  MHz.  The  general 
ciualitative  features  described  in  Figure  8  are  seen  to  be  present 
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FIGURE  8.  SCHEMATIC  DIAGRAM  OF  A  TYPICAL  EXPERIMENTAL  CURVE 
FOR  THE  REFLECTED  AMPLITUDE  R(6)  AND  PHASE  <j>(6) 

OF  THE  REFLECTED  WAVE. 


FIGURE  9.  EXPERIMENTAL  ARRANGEMENT  OF  THE  SOURCE  AND  RECEIVER. 


FIGURE  10.  EXPERIMENTAL  MEASUREMENTS  ON  A  I  T  TAN  I  I'M  ALLOY 

AT  i.l  MHZ  ( 2  5  °c: )  SHOVING  THE  REFLECTED  AMPLITUDE 
(RELATIVE  UNITS)  K(-~)  IN  THE  UPPER  CURVE  AND  THE 
PHASI  :('•)  LOWER  CURVE.  IN'  THESE  CURVES,  THF.  TIC 
MARKS  ARE  r  AND  THE  RANGE  OF  THE  PHASE  DATA  IS 
3HOr.  Till:  DATA  AT  THE  RIGHT  SUPPLIES  FURTHER 
INFORMATION  REGARDING  THE  SPECIMEN. 
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FIGURE  II.  SAME  SPECIMEN  AS  SHOWN  IN  FIGURE  10,  BUT  WITH  A 

BETTER  MINIMIT1  IN  R((>).  DISTORTIONS  IN  THE  PHASE 
CURVE  .;(•')  (BOTTOM  CURVE)  ARE  READ  I EY  APPARENT. 
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in  this  measured  data. 


However,  a  considerable  number  of  "anomalies"  are  also 
encountered  in  making  these  measurements.  When  we  are  careful  to 
locate  the  frequency  of  least  reflection  (FLR)  and  when  we  have  very 
carefully  oriented  the  sample  (see  Experimental  Procedure  Section  4.3), 
we  often  observe  a  "distorted"  phase  curve  for  the  same  sample 
(see  Figure  11) . 

In  the  earlier  experiments,  we  were  inclined  to  attribute  these 
phase  curve  "distortions"  to  experimental  error  since  they  were 
evidently  very  sensitive  to  set  up  parameters.  However,  after  we 
introduced  temperature  controls  and  a  gimbaled  sample  holder,  it  was 
shown  that  results  were  highly  reproducible.  That  is,  a  sample  could 
be  measured,  removed,  and  then  remeasured  with  nearly  identical 
results.  Thus,  we  found  that  at  the  same  frequency,  temperature, 
sample  elevation,  and  orientation  (and  direction  of  examination  along 
the  surface)  we  obtained  essentially  identical  amplitude  curves. 

Slight  changes  in  the  orientation  of  the  normal  to  the  sample 
surface  could  result  in  the  appearance  of  distinct  distortions  in 
the  phase  curves  at  the  minimum.  We  soon  noticed  that  when  these 
distortions  occurred,  we  invariably  had  produced  the  least  minimum 
in  the  amplitude  curves.  That  meant  that  we  had  in  fact  an  even 
better  set  up  than  the  original  measurement  (using  standard  experimental 
procedure  of  setting  up  on  45°  and  maximizing  the  reflected  amplitude 
by  adjusting  the  sample  orientation)  which  had  produced  an  "undistorted" 
phase  curve. 
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At  this  point  in  our  experiments,  we  decided  to  examine  the 
actual  reflected  signals  at  the  minimum  amplitude  location.  If 
this  signal  contains  some  unwanted  component,  it  could  explain  the 
anomalous  phase  curves.  In  Figure  1 2  we  illustrate  in  schematic 
form  the  typical  waveforms  observed  near  the  minimum  of  the  reflected 
amplitude. 

Signals  of  this  type  are  invariably  seen  in  numerous  samples 
whenever  a  distorted  phase  curve  is  seen.  By  slightly  changing  the 
sample  orientation  until  the  signal  amplitude  increases,  the  signal 
distortions  become  less  obvious  to  the  eye  and  a  "good"  phase  curve 
is  again  obtained.  This  explains  why  the  sample  orientation  is  so 
important  in  observing  this  effect. 

There  is  clearly  a  second  harmonic  contamination  present  in 
these  signals.  We  have  observed  that  as  the  primary  reference 
frequency  is  increased,  the  harmonic  frequency  increases  accordingly. 
Because  this  second  harmonic  contamination  is  of  comparable  amplitude 
to  the  incident  wave  (at  9C)  at  the  minimum  of  R(9) ,  it  was  at  first 
thought  that  this  was  the  cause  of  the  distorted  phase  curves. 

However,  the  vector  volt  meter  (phase  and  amplitude  measuring 
system)  we  are  employing  was  designed  to  be  insensitive  to  a  second 
harmonic  contamination  of  this  magnitude  (it  is  filtered  out).  Thus, 
in  spite  of  the  very  obvious  second  harmonic  contamination,  it  is 
not  responsible  for  the  observed  phase  curve  distortions. 

By  band-pass  filtering  the  received  signal,  it  was  eventually 
found  that  these  phase  distortions  are  due  to  signals  at  the  original 


Time 


SCHEMATIC  DIAGRAM  SHOWING  THE  OBSERVED  SIGNAL  AT  THE 
MINIMUM  OF  R(0) .  THIS  SIGNAL  APPEARS  TO  BE  THE  SUM 
OF  A  PRIMARY  WAVE  OF  FREQUENCY  w0  AND  A  SECOND  HARMONIC 
FREQUENCY  2ojq  AT  90°  TO  THE  PRIMARY.  WHENEVER  SUCH  A 
SIGNAL  WAS  OBSERVED,  THE  PHASE  CURVES  WERE  INVARIABLY 
DISTORTED.  HOWEVER,  IT  WAS  SOON  LEARNED  THAT  THIS 
SECOND  HARMONIC  CONTAMINATION  WAS  NOT  RESPONSIBLE  FOR 
THESE  DISTORTIONS  SINCE  THE  PHASE  METER  DID  NOT 
RESPOND  TO  THIS  COMPONENT  (SEE  TEXT  FOR  A  DISCUSSION). 
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reference  frequency  u>q. 

We  tentatively  conclude  that  this  additional  component  at 

the  reference  frequency  (jl)q  is  the  result  of  some  effect  due  to 

transmission  through  the  water,  a  sample  dependent  effect,  or  both. 

We  note  that  nonlinear  interactions  could  be  responsible  either  in 

30-34 

the  water  or  the  sample  or  both.  If  this  is  the  case,  the  signal 

at  the  reference  frequency  ojq  would  have  to  interact,  say  with  a 
second  harmonic  at  2u)  ,  to  yield  a  signal  at  2oj  -  to  =  u  .  It  is 
very  difficult  to  say  at  this  time  just  what  the  mechanism  is. 

Other  possibilities  still  cannot  be  eliminated.  We  are  certain, 
however,  that  these  effects  are  not  due  to  electronic  contamination 
or  other  peculiarities  of  our  experimental  apparatus.  These  effects 
are  definitely  occurring  along  the  sound  path  and  there  are  indications 
that  it  is  the  sample  and  not  the  water  path  that  is  introducing 
this  new  signal.  It  also  begins  to  look  as  if  this  phase  distorting 
signal  only  appears  when  we  are  at  a  critical  angle. 

To  see  this,  we  will  now  discuss  the  second  kind  of  anomaly 
that  is  frequently  observed  on  the  same  sample.  In  Figure  11  it  may 
be  noted  that  just  to  the  left  of  the  minimum  in  R(9)  (about  2  1/2° 
or  1/4  tick  mark  on  the  photo),  there  is  a  slight  indication  of  a 
second  minimum  since  the  average  slope  is  smaller  in  this  region. 

We  discover  in  fact  that  by  rotating  the  specimen  90°  (smaller 
angles  also  work)  in  the  plane  of  its  surface  about  its  normal 
vector  we  obtain  the  measurements  shown  in  Figure  13(a). 

In  this  figure  there  are  clearly  two  dips  and  two  associated 
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"distorted"  phase  shifts.  Again  it  should  be  stressed  that  measurements 
of  this  sort  are  repeatable,  nor  are  these  multiple  dips  unusual.  We 
see  indications  on  many  samples  of  a  second  dip.  Two  of  the  most 
dramatic  results  being  those  found  in  CaF2  and  a  4340  steel  (see 
Figures  13(b)  and  (c)).  It  is  clear  from  Figure  13  that  the  mechanism 
causing  the  phase  distortions  (whatever  its  exact  nature)  is  present 
at  each  dip  but  does  not  seem  to  be  present  elsewhere.  This  strongly 
implies  that  the  effect  (phase  distortions)  is  sample  dependent. 

Apart  from  these  observations  regarding  phase  distortions, 
the  very  existence  of  multiple  dips  is  a  surprising  and  completely 
unexpected  phenomenon.  None  of  the  existing  theories  predict  such 
an  effect.  A  preliminary  interpretation,  which  is  by  no  means  the 
only  possible  interpretation,  is  that  a  surface  wave  is  being 
generated  at  each  of  these  amplitude  dips.  If  this  were  the  case, 
we  must  definitely  reexamine  the  theoretical  foundations  since 
conventional  theories  with  conventional  boundary  conditions  predict 
only  one  surface  wave  along  a  liquid-solid  interface  (isotropic 

3  S 

solid) . 

In  Table  I,  we  list  the  results  for  8^  at  the  frequency  of 


least  reflection  on  several  different  types  of  materials. 
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TABLE  I:  PRINCIPAL  CRITICAL  ANGLE  0C  FOR  VARIOUS  MATERIALS  AT  THE 
FREQUENCY  OF  LEAST  REFLECT  FLR 


Metals  and  Alloys 


M°) 

FLR (MHz) 

Copper 

47.5 

4.46 

Brass 

50.7 

2.14 

Aluminum(707  5) 

33.5 

4.29 

Titanium 

31.6 

3.10 

4340  Steel 

33.8 

2.50 

r 

Woods  Metal 

— 

-  (no  0C  15-55°) 

"Ceramics" 

Tungsten  Carbide 

- - 

-  (no  0C  15-55°) 

Tungsten  Powder  Composite 

37.8 

4.9 

Alumina (Coors) 

17.8 

4.9 

Polymers 


Rubber  -  - (no  0C  15-55°) 

Polys tryene  -  - 


Glasses 


Obsidian  26.9  2.62 

Silica  26.5  2.55 
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Complex  Silicates 


ec(°) 

FLR(MHz) 

Granite(fine  grained) 

26.2 

4.00 

"Tiger  Eye" 

25.2 

4.00 

Epidote 

25.0 

5.00 

Single  Crystals 

Silicon  (Si) 

21.1 

4.45 

Saphire(Al203) 

15.6 

3.79 

Fluorite(CaF2) 

28.5 

4.77 

QuartE(Si02) 

28.6 

3.74 

4 

4 
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4 . 6  Frequency  of  Least  Reflection  (FLR) 


An  examination  of  the  data  of  Table  I  shows  that  the  measured 
frequencies  of  least  reflection  are  typically  less  than  5  KHz  for  a 
wide  range  of  material  types. 

9 

This  observation  was  somewhat  surprising  since  other  authors 
report  FLR's  typically  above  10  MHz.  Richardson  has  also  indicated 
(private  communication)  that  in  experiments  performed  by  Becker  and 
Richardson,  the  minimum  at  FLR  did  not  appear  to  be  as  deep  as  we 
generally  find.  We  have  not  yet  determined  the  cause  of  these 
differences  but  since  our  electronics  responds  well  in  the  range 
where  measurements  were  made,  we  are  confident  that  these  FLR's  are 
correct.  It  is  possible  that  some  aspect  of  our  experimental 
apparatus,  such  as  the  focused  lens,  is  responsible  for  these 
differences.  However,  it  does  not  seem  likely  that  a  function  which 
undergoes  a  minimum  (R ( 0 ) )  and  which  depends  on  material  properties 
will  be  effected  by  a  different  measuring  procedure.  Ihe  rrinimum 
will  still  occur  at  a  frequency  oj(FLR)  which  should  be  essentially 
independent  of  the  apparatus  used  to  observe  it. 

For  these  reasons,  our  experimentally  determined  FLR's  seem 
to  be  at  odds  with  those  found  by  other  investigators  although  we 
have  not  examined  the  same  samples  as  these  investigators. 
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').0  PROSPECTS  FOR  COMPARISON  OF  EXPERIMENTS  WITH  EXISTING  THEORIES 

A  number  of  authors  working  independently  have  arrived  at 
rather  similar  theoretical  descriptions  for  the  critical  angle 
phenomenon  in  spite  of  what  at  first  appears  to  be  different 
assumptions . 

l  j 

We  will  discuss  here  the  work  of  Richardson  and  that  of 

24 

Pitts  as  representative  of  the  available  theories.  Various 
qualitative  predictions  made  by  these  models  will  be  discussed  and 
characterized.  Comparison  between  these  predictions  and  our 
experimental  results  will  be  made.  Questions  will  then  be  raised 
regarding  possible  modifications  of  these  theories. 

5 . 1  The  Richardson  Model 

It  has  been  known  experimentally  for  some  time  that  at  a 

critical  angle  both  an  amplitude  dip  and  a  phase  shift  occur  in  a 

reflected  acoustic  (or  elastic)  wave  incident  on  a  boundary  between 

6-8 

two  solids  (or  a  solid-liquid  boundary).  The  mechanism  for  this 

behavior  is  the  conversion  of  the  incident  longitudinal  or  pressure 
wave  into  a  surface  or  Rayleigh-type  wave  that  propagates  along 
the  boundary  and  extends  roughly  one  wavelength  into  the  media  on 
either  side  of  the  boundary.  This  conversion  of  incident  wave 
energy  to  surface  wave  energy  and  the  subsequent  reradiation  back 
into  the  water  accounts  for  the  dip  in  the  amplitude  of  the 


reflected  wave. 


Until  recently,  however,  the  theory  for  this  behavior  was 
unsatisfactory.  Whereas  surface  wave  solutions  have  long  been  known 

3  5 

to  exist  in  classical  "hookean"  materials,  the  mechanism  for 
exciting  these  waves  was  not  well  understood  until  attenuation  was 
added  to  the  analyses.  The  Richardson  linear  visco-elastic  wave 
model  among  others  has  updated  the  classical  Hooke's  law  theory  by 
including  the  effects  of  attentuation  via  the  addition  of  "viscous" 
terms  in  the  stress-strain  relation.  It  has  been  found  that  these 
viscous  terms  indirectly  determine  the  coupling  between  an  incident 
wave  and  the  respondent  surface  wave.  Without  such  terms,  very  little 
if  any  surface  wave  energy  will  result  at  the  boundary  between  two 
media. 

Richardson's  linear  visco-elastic  model  in  a  one-dimensional 
approximation  may  be  described  by  the  stress-strain  relation 

a  =  KiC  +  —  e  (6) 

U) 

where  higher  order  terms  involving  £  etc.  are  assumed  to  be  unimportant, 
and  where  the  "viscosity"  coefficient  r)  =  K2/U)  is  frequency  dependent. 
This  is  another  way  of  saying  that  the  "response  time”  x  is  distributed 
as  1/u) 


T 


k2 

K  ho 


(7) 


Using  this  simple  model,  it  is  easily  shown  that  the  stress-strain 


ellipse  is  frequency  independent  since  we  find  for  £  =  sin  wt 


--  36 

Frequency  independence  of  the  a  -  £  ellipse  was  shown  by  Knopoff 

to  be  an  observed  property  of  many  materials  over  a  very  wide  range  of 

1 1 

acoustic  frequencies.  Richardson  and  Becker  have  shown  experimentally 
that  this  frequency  independent  behavior  also  holds  for  metals  and 
other  similar  materials  in  the  MHz  range  of  frequencies. 

By  generalizing  Equation  (4)  to  a  tensor  equation  in  three 
dimensions 


,  .  Dkl 

„kl  ,  lj 

ij  kl  to  kl,t 


Richardson  was  able  to  derive  the  modified  elastic  wave  equations 

describing  the  propagation  of  P  and  S  waves  in  isotropic,  homogeneous 

media.  From  these  equations,  he  was  then  able  to  solve  the  particular 

problem  of  the  reflection  of  an  incident  wave  from  the  boundary 

between  two  dissimilar  media.  In  particular,  he  found  the  reflection 

1  3 

coefficient  at  a  liquid-solid  boundary.  Working  independently, 

1  5 

Merkulova  arrived  at  the  same  result.  The  essence  of  the  theory 
resides  in  the  equation  for  the  reflection  coefficient  found  by  these 
authors 


cosz2y22  +  (C22/C12)  sin2Yi2sin2Y22  +  p i c 1 1 siny 1 2 /p2c 1 2sinyi t 


cos  2y 2 2  +  (c22/ci2)zsin2Yi2sin2Y22  -  P jc 1 j sinYi 2/ P2C 1 2sinYi i 


where  C12  =  velocity  of  longitudinal  waves  In  the  solid. 


C22  =  velocity  of  shear  waves  in  the  solid, 

cii  =  velocity  of  longitudinal  waves  in  the  liquid, 

0  =  Y11  =  angle  of  incidence  of  the  longitudinal  wave  in  the  liquid, 

Y12  =  angle  of  refracted  longitudinal  wave  in  the  solid, 

Y22  =  angle  of  refracted  shear  wave  in  the  solid. 

Pi  =  density  of  liquid 
pz  =  density  of  solid. 

Angles  are  measured  in  the  counterclockwise  direction  from  the  x-axis 
(plane  of  the  interface).  The  velocities  are  the  complex  quantities 

c  =  v  ( 1  +  iva/co),  (11) 

v  =  normal  velocity, 
a  =  attenuation, 
w  =  radian  frequency. 

The  refracted  angles  are  found  by  the  generalized  Snell's  law, 


cosy  •  . 

— =  h  .  (12) 

c  .  . 
ij 

These  equations  have  been  specialized  to  a  two  layered 
liquid-solid  medium  although  the  theory  holds  for  solid-solid 
interfaces  as  well.  This  model  is  said  to  predict  observed  behavior 
for  unbounded  beams  with  good  accuracy.  Both  the  dip  in  amplitude 
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REFLECTED  AMPLITUDE  VERSUS  ANGLE  OF  INCIDENCE  FOR  A  WATER/ 
STAINLESS  STEEL  INTERFACE.  THE  OBSERVED  REFLECTION  FACTOR 
(SOLID  LINE)  DIFFERS  FROM  THAT  CALCULATED  BY  NONATTENUATION 
WAVE  MODELS  (DASHED  LINE)  IN  THE  AREA  OF  THE  RAYLEIGH  CRIT¬ 
ICAL  ANGLE. 
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COORDINATE  SYSTEM 
FOR  RICHARDSON'S  CALCULATION 
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CALCULATED  AND  MEASURED  PHASE  AND  AMPLITUDE  OF  REFLECTED 
WAVES  AT  A  WATER- 304  STAINLESS-STEEL  INTERFACE  FOR  FRE¬ 
QUENCIES  OF  10.0,  14.5,  AND  18.0  MHz  FOR  RICHARDSON'S  MODEL. 
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and  the  phase  change  are  predicted  quantitatively  as  well  as 
qualitatively . 

In  addition  to  describing  the  behavior  of  the  reflected 
wave  at  the  critical  angle,  the  theory  indicates  that  there  exists  a 
frequency  at  which  there  is  zero  reflection  at  the  critical  angle 
called  the  frequency  of  least  reflection  (FLR).  This  point  is  said 
to  be  highly  dependent  upon  grain  size,  and  appears  to  present  a 
method  for  its  measurement. 

The  theory  developed  by  Richardson  is  nearly  complete  with 

respect  to  the  steady  state  solution.  In  the  course  of  the  experimental 

verification,  however,  some  interesting  transient  effects  were 
1  ? 

observed.  It  was  noted  that  the  phase  response  near  the  critical 
angle  involved  a  time  lag.  That  is,  it  took  a  finite  amount  of  time 
after  initiation  of  the  ultrasonic  burst  before  the  phase  of  the 
reflected  wave  changed  from  the  classically  predicted  value  to  that 
obtained  with  the  Richardson  theory.  Thus,  it  appears  that  extension 
of  the  theory  to  include  transient  effects  might  be  profitable. 

5 . 2  The  Model  of  Pitts 

Whereas  Richardson  has  provided  a  particular  model  based  on  a 
fundamental  stress-strain  relation  including  attenuation,  Pitts  at 
first  appears  to  take  a  different  and  more  conventional  view  by  intro¬ 
ducing  attenuation  directly  in  the  wave  numbers.  The  end  results  of 
these  two  approaches  are  nevertheless  rather  similar  and  may  ultimately 
involve  identical  physical  models  in  some  cases.  We  will  review 
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rather  closely  the  derivation  given  by  Pitts  in  his  1976  thesis 


and  discuss  the  approach  as  we  proceed. 

Pitts  remarks  that  his  derivation  follows  the  standard 
approach  given  in  most  textbooks  although  his  results  differ  slightly 
from  those  of  Brekhovskikh  for  example.  In  what  follows,  we  will 
employ  a  set  of  symbols  that  are  different  than  those  employed  by 
Richardson.  Richardson's  formulation  does  not  appear  to  use 
"conventional"  notation.  For  example,  he  employs  complex  velocities 
whereas  Pitts  and  others  employ  complex  wavenumbers.  These  differences 
may  make  it  difficult  to  compare  the  results  of  these  authors,  but 
such  differences  alone  do  not  lead  to  formally  different  answers. 


waves  (s=shear,  d=longitud inal  in  solid).  The  law  of  reflection  is 


assumed  to  hold,  and  the  propagation  directions  of  the  sound  waves  are 
assumed  to  be  given  by  Snell's  law 


Here,  V  is  the  longitudinal  wave  velocity  in  water,  is  the  long¬ 
itudinal  wave  velocity  in  the  solid,  and  Vg  is  the  shear  wave  velocity. 

Although  attenuation  is  not  explicitly  introduced  at  this 
stage,  these  velocities  and  angles  should  be  assumed  to  be  generally 
complex  to  account  for  the  attenuation. 

Pitts  then  writes  down  the  particle  displacement  vectors 


U 


-  ik*r  -iojt 
Ae  e 


(14) 


for  the  two  media. 

For  medium  1  (the  liquid): 

Aiexp[ik(xsin0  +  zcos0)-iu)t] 

+  B iexp [ ik(xsin0  -  zcos0)-iwt] 

For  medium  2  (the  solid): 

U2  =  A  , exp [ ik,  (xsin9  ,  +  zcosB  , )-iu>t  ] 
dad  d 

+  B ,exp[ik , (xsin0 .  -  zcos0.)-iwt] 
_d  d  d  d 

+  A  exp[ik  (xsin9  +  zcos0  )-iiotl 
s  s  s  s  4 

+  B  expfik  (xsinG  -  zcosQ  )  —  lout 7 
s  s  s  s 

The  sound  wave  amplitudes  are: 

Ai  =  Ai[cos0k  +  sinOij 
Bi  =  Bi[-cos0k  +  sinBi] 

A,  =  A,[cos0,k  +  sin0.il 

d  d  d  d 

B.  =  B,  f-cosO.k  +  sin0,i] 

d  d  d  d 

A  =  A  [sin0  k  -  cos0  i] 
s  s  S/>  s^ 

B  =  B  [sinO  k  +  cosO  i] 
s  s  s  s 


(15) 


(16) 


(17) 
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where  the  unit  direction  vectors  x  =  i  and  z  =  k.  The  symbol  k 

A  A 

should  not  be  confused  with  a  wavevector  k  and  i  should  not  be 
confused  with  i  =  \J- 1 . 

The  following  conventional  boundary  conditions  on  the  stresses 
and  particle  velocities  are  then  applied  to  these  solutions. 

_  /\ 

A)  The  normal  component  of  the  particle  displacement  (U*k) 
is  continuous  across  the  boundary 

B)  The  normal  component  of  stress  is  continuous  across  the 
boundary 

C)  The  shear  stress  is  zero  at  the  boundary 

After  considerable  manipulation  (not  included  in  Pitts' 
thesis),  he  finds  the  following  reflection  coefficient  for  the  unbounded 
plane  wave 

\s(kx)  =  N/D  (18) 


where  k^  is  the  x-component  wavenumber  of  the  incident  wave  and 
where 


and 


N  =  (k2  -  2k2)2  +  4k2K  K,  -  pkMKd/K 
s  x  x  s  d  s 

D  =  (k2  -  ?k2 ) 2  +  4k2K  K,  +  pk4K,/K 
s  x  x  s  d  s  d 


p  =  p  liquid/p  solid 

k  =  (u)/v)  (1  +  i  A/ 2ti  ) 

k  =  (w/v  )(1  +  iA  /  2tt ) 
s  s  s 

kd  =  (w/vd)(1  +  iAd/2TT) 
K  =  (k2  -  k2)l/2 

K .  =  (k2  -  k2)1/2 
d  d  x 

K  =  (U7  -  k 2  ">  1  / ' 


(19) 


(20) 
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Note  that  the  loss  parameters  A,  A^ ,  Ag  appear  in  the  wave- 

numbers  k,  k , ,  k  rather  than  in  the  velocities  V,  V  ,  V„.  Therefore 
d  s  s  d 

v,  v  ,  v.  are  real  and  V,  V  ,  V,  are  complex, 
s  d  s  d 

The  function  R  (k  )  depends  both  on  the  incident  angle  0 

Lj  X 

and  the  frequency.  The  frequency  dependence  enters  because  A,  A  ,  A, 

s  d 

are  generally  frequency  dependent. 


5 . 3  Comparison  of  Pitts'  and  Richardson's  Models 

Although  the  symbols  employed  in  Pitts  formulation  of  R(k^) 
(Equation  18)  are  different  from  those  employed  by  Richardson 
(Equation  10),  these  expressions  lead  to  very  similar  results  and  are 
based  on  rather  similar  ideas  even  though  they  may  seem  to  represent 
different  approaches. 

1)  Both  Pitts  and  Richardson  employ  the  same  boundary 
conditions. 

2)  Although  Richardson  introduces  attenuation  at  the  level 
of  Hooke's  law  and  then  seeks  solutions  to  a  modified 
equation  of  motion,  his  result  should  not  be  formally 
different  from  a  derivation  such  as  Pitts  where  solutions 
which  include  attenuation  are  postulated  ad-hoc.  Both 
approaches  would  seem  to  be  satisfactory  and  could  very 
well  lead  to  identical  results.  Both  approaches  are, 
however,  strictly  phenomenological. 

3)  A  single  dip  in  the  amplitude  and  an  associated  phase  shift 

are  predicted  by  both  models  at  the  critical  angle.  Both 

models  predict  that  the  phase  shift  is  insensitive  to 

changes  in  shear  wave  attenuation  whereas  the  reflection 

factor  (R(k  ))  is  predicted  to  be  very  sensitive  to  the 

choice  of  A  . 

s 

4)  Both  models  predict  the  phenomenon  known  as  the  frequency 
of  least  reflection  (FLR). 

To  summarize,  we  can  say  that  the  Richardson  model13  which 


I 


1  5 

produces  results  identical  to  those  of  Merkulova  also  produces 

2  4.  2  5  2  6-2  9 

results  qualitatively  very  similar  to  those  of  Pitts  and  others  who 

have  investigated  similar  models.  Although  detailed  comparisons  have 
yet  to  be  made,  there  would  seem  to  be  little  if  any  essential  physical 
differences  between  these  models  in  so  far  as  predictions  are  concerned. 

From  our  experimental  results,  there  are  indications  that  all 
of  these  models  may  have  drawbacks.  In  this  context,  the  differences 
between  one  model  and  a  very  similar  model  are  less  significant  than 
the  apparent  inability  of  any  model  to  explain  experimental  data. 

5. A  Theoretical  Difficulties 

The  analysis  of  Richardson,  Pitts,  and  other  authors  that  have 
investigated  this  problem  employ  conventional  boundary  conditions  i.e., 

A)  Normal  component  of  particle  displacement  is  continuous 
across  the  boundary  (liquid-solid) 

B)  Normal  component  of  stress  is  continuous  across  the  boundary 

C)  The  shear  stress  is  zero  at  the  boundary. 

It  is  also  tacitly  assumed  by  these  authors  that  one  can 
examine  the  ideal  case  (Hooke's  law  behavior)  and  by  appropriate 
modifications,  such  as  the  addition  of  attenuation,  (either  in  a 
modified  Hooke's  law  response  or  in  the  wavenumbers)  obtain  the  actual 
response  of  a  material  or  boundary.  This  of  course  is  the  traditional 
phenomenological  approach  followed  in  most  treatises  on  the  theory 
of  elasticity.  However,  this  traditional  approach  raises  a  number  of 
serious  questions  of  a  physical  nature. 
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Consider  the  standard  problem  of  computing  the  reflection  and 


transmission  coefficients  from  a  plane  boundary  between  two  infinite 
elastic  media.  In  particular,  consider  the  case  of  a  liquid-solid 
interface. 

Brekhovskikh  ^  like  many  other  authors  in  the  past,  has  derived 
these  coefficients  for  the  lossless  case  and  demonstrated  that  the 
energy  of  the  elastic  waves  is  conserved.  These  authors,  however, 
have  not  apparently  gone  on  to  ask  if  linear  momentum  is  conserved. 
Consider  the  reflection  and  transmission  coefficients  as  derived  by 
Brehhovskikh 1  which  is  illustrated  graphically  in  his  book  on  page  46 
and  reproduced  in  Figure  14. 

While  it  may  be  satisfying  that  energy  is  conserved  in  this 
process,  a  serious  question  can  be  raised  regarding  the  correctness 
of  this  analysis  on  the  basis  of  momentum  conservation. 

Note  that  both  the  transmitted  shear  wave  and  the  transmitted 
longitudinal  wave  vanish  at  the  longitudinal  critical  angle  0^.  At 
this  angle,  there  is  100%  reflection  of  the  incident  wave  (E'/E  =  1). 
Similarly,  at  the  shear  critical  angle  8g ,  there  is  100%  reflection 
of  the  incident  wave  E'/E  =  1  (see  Figure  14).  Beyond  0g,  of  course, 
we  have  the  additional  possibility  of  the  Rayleigh  critical  angle 
phenomenon. 

Consider  then  a  wave  which  is  100%  reflected  from  a  solid  as 
shown  in  Figure  15.  Let  us  calculate  the  change  in  momentum  of  this 
incident  sound  wave.  The  x-components  of  momentum  will  clearly  be 
unaltered;  however,  there  will  be  a  finite  change  in  the  z-component. 
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FIGURE  14.  THE  CONVENTIONAL  ANALYSIS  OF  THE  REFLECTION  AND 

TRANSMISSION  COEFFICIENTS  OF  A  LIQUID- SOLID  BOUNDARY 
UP  TO  THE  SHEAR  CRITICAL  ANGLE  8S .  CURVES  TAKEN 
FROM  BREKHOVSKIKH.1  0L  IS  THE  LONGITUDINAL  CRITICAL 
ANGLE.  IN  THE  LOSSLESS  CASE,  ENERGY  IS  CONSERVED  AND 


FIGURE  15.  A  TOTALLY  REFLECTED  SOUND  WAVE.  HOW  IS  THE  MOMENTUM 
TRANSFERRED  TO  THE  SOLID  EVEN  IF  ENERGY  IS  CONSERVED? 


If  the  wave  is  essentially  unbounded,  the  group  and  phase 
velocities  are  the  same.  If  v  is  the  phase  velocity  in  the  liquid, 
the  z-component  of  momentum  of  the  incident  sound  wave  will  be 
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proportional  to 

P  (i)  «  -  —  cosB  .  (21) 

z  v 


Similarly,  the  z-component  of  the  reflected  wave  will  be  proportional  to 


.\P  a  COSO  . 

Z  V 

Clearly  since  0  is  generally  an  angle  different  from  tt/2,  i.e. 

=  6  or  G^,  this  quantity  is  not  zero. 

By  the  law  of  conservation  of  linear  momentum,  the  change  in 
momentum  of  the  solid  must  be  proportional  to  -AP.,.  Therefore, 
momentum  must  be  transferred  to  the  solid  from  the  incident  wave. 

This  is  clearly  impossible . since  no  bulk  waves  exist  in  the  solid  at 
the  angles  0  =  6g  or  0^. 

Thus,  the  conventional  analysis  of  the  lossless  case,  using 
conventional  boundary  conditons,  leads  to  a  fundamental  error  and 
casts  doubt  on  other  similar  analyses  including  the  aforementioned 
analyses  of  the  critical  angle  phenomenon. 

5 . 5  Some  Possible  Remedies 

To  illustrate  how  serious  the  foregoing  paradox  is,  one  need 
only  recall  that  the  same  analysis  carried  out  for  a  totally  reflected 
light  beam  shows  that  the  momentum  transferred  upon  reflection 

,37 

varies  as  E/e  .  Hence,  a  sound  beam  having  the  same  energy  density 

as  a  light  beam  carries  c/v  times  as  much  momentum.  It  is  acceptable 

in  most  optical  applications  to  ignore  the  small  momentum  transfers 
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upon  reflection  (except  for  very  intense  beams  )  but  one 
clearly  cannot  ignore  the  relatively  large  momentum  transfers  involved 
in  acoustic  reflections.  How  then  is  the  traditional  analysis  of 
Brekhovskikh  and  others  to  be  modified  in  order  to  correct  such 
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deficiencies?  This  question  has  not  as  yet  been  answered,  however, 
we  have  been  able  to  characterize  the  problem  sufficiently  to 
identify  the  probable  source  of  these  conceptual  difficulties. 

First,  the  practice  of  obtaining  approximate  solutions  by 
examining  a  lossless  case,  as  Brekhovskikh  does,  has  a  very  questionable 
validity.  We  note  that  sound  waves  are,  strictly  speaking,  impossible 
unless  nonlinearities  (interactions)  and  their  inevitable  associated 
losses  are  present  (see  Appendix  I).  This  is  not  a  trivial  point 
since  it  means  that  even  in  principle  there  is  no  lossless  case  in 
acoustics.  Such  a  situation  does  not  even  exist  in  a  first  approxi¬ 
mation.  Nor  is  it  possible,  strictly  speaking,  to  describe  these 
processes  with  linear  partial  differential  equations. 

Approaches  based  on  these  oversimplified  phenomenological 
models  (lossless,  linear  systems)  may  give  reasonable  agreement  with 
experience  in  some  situations,  but  we  should  not  be  too  surprised 
if  they  fail  in  many  others. 

A  real  sound  wave  carrying  energy  and  momentum  incident  on  a 
solid,  say  at  a  liquid-solid  interface,  must  eventually  transfer  an 
amount  of  energy  and  momentum  to  the  so  called  ’'zero  frequency" 
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modes  of  the  solid.  Such  zero  frequency  modes  represent  the  motion 
of  the  solid  as  a  whole.  A  solid,  immersed  in  a  liquid  in  a  state  of 
free-fall  and  subject  to  an  incident  sound  pulse  (duration  t) , 
would  thus  begin  to  translate  with  a  constant  average  velocity  V 
(neglecting  resistance  of  the  fluid)  just  after  the  pulse  had  left 
the  object.  Clearly,  there  will  be  an  acceleration  of  the  body  of 
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order  V/i.  The  process  of  transferring  acoustic  momentum  to  the 
solid  is  not  a  linear  one  since  a  wave  incident  at  an  average 
frequency  of  u>  must  eventually  he  converted  to  these  zero  frequency 
mod  es . 

In  the  language  of  solid  state  physics,  we  can  say  that  the 
initial  phonon  spectrum  of  the  solid  body  has  been  converted  to 
another  spectrum  raising  the  temperature  of  the  solid  object.  In 
tile  process,  the  zero  frequency  modes  have  been  excited. 

What  the  foregoing  discussion  suggests  is  that  the  traditional 
b o u nd a rv  conditions  employed  in  such  problems  must  be  altered  in 
order  to  obtain  a  consistent  theory.  The  reason  for  this  conclusion 
is  the  fact  that  the  interaction  at  a  boundary  (liquid-solid)  is  a 
local  one.  At  the  instant  a  real  sound  wave  impinges  on  the  boundary, 
it  begins  the  process  of  transferring  momentum.  This  can  take  place 
in  two  principal  ways  on  a  macroscopic  scale: 

A)  Via  wave-like  motions  (wave  equations) 

B)  Diffusion-like  motions  (damping,  viscosity)  diffusion  equations 

C)  Combinations  of  A)  and  B) . 

Wave-like  motions,  of  course,  refers  to  elastic  waves,  and 

diffusion-like  motions  refers  to  various  viscosity  terms  and  leads  to 
1  ? 

transient  effects.  Now  as  mentioned,  most  authors  describing  the 
critical  angle  phenomenon  have  correctly  recognized  that  the  critical 
angle  phenomenon  depends  in  some  way  on  the  presence  of  attenuation. 
Without  attenuation,  the  phenomenon  is  certainly  inexplicable.  These 
explanations  cannot  be  complete,  however,  since  they  employ  traditional 
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boundary  conditions.  There  is  no  indication  that  they  have  resolved 
the  momentum  problem  for  example. 

Any  narrow  region  where  the  solution  of  a  differential  equation 
changes  rapidly  is  referred  to  in  advanced  mathematical  treatises 
as  a  boundary  layer.  It  may  well  be  that  the  complete  solution  we 

h  4 

are  seeking  will  involve  solutions  based  on  Boundary-Layer  Theory. 

To  summarize,  we  can  say  that  a  correct  treatment  of  the  Rayleigh 
critical  angle  phenomenon  will  also  automatically  provide  a  correct 
general  expression  for  scattering  at  all  angles;  not  just  the  Rayleigh 
critical  angle.  In  particular,  the  model  will  be  consistent  with 
energy  and  momentum  conservation  in  the  presence  of  nonlinearities 
and  losses.  Such  a  model  must  also  explain  the  phenomenon  of  multiple 
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6.0  BRIEF  REVIEW  OF  PHASE  II  PLANS 


The  Phase  I  study  has  delineated  a  number  of  problems  with 
the  conventional  phenomenological  descriptions  of  the  so  called 
Rayleigh  critical  angle  phenomenon.  Accordingly,  the  Phase  II 
program  in  brief  outline  will  consist  of  the  following  elements: 

Task  1:  Theoretical  Development 

Continued  Studies  in  Basic  Elasticity  Theory 
Damping 

Nonlinearities 
Second  Order  Effects 
Anisotropy  and  Discontinuities 

Coupling  Mechanisms 

Momentum  Balance  (see  Figure  16 ) 

Study  of  Focused  Beams 

Other  Focal  Region  Effects 

Electromagnetic  Properties 
Induced  Temperature  Gradients 
Applied  Stress  Effects 

Task  2:  Experimental  Backup  for  Task  1 
Materials  with  Double  Dips 
Monitor  Focal  Region  (see  Figure  17) 

Experiments  in  Momentum  Balance 


Higher  Frequency  Experiments 


z 


Incident  Wave 


Reflected  and  Reradiated 
Components 


FIGURE  16.  THIS  FIGURE  ILLUSTRATES  THE  CONVENTIONAL  PARTITIONING  OF 
AN  INCIDENT  WAVE  INTO  A  RESPONDENT  "LEAKY"  SURFACE  WAVE 
AT  THE  CRITICAL  ANGLE  AND  REFLECTED  AND  RERADIATED 
COMPONENTS.  ONE  CAN  EASILY  IMAGINE  THAT  MOMENTUM  IS 
CONSERVED  ALONG  X  BUT  UNLESS  THERE  IS  A  BULK  WAVE  IN  THE 
SOLID,  HOW  CAN  IT  BE  CONSERVED  ALONG  Z?  IF  THERE  ARE  NO 
BULK  WAVES  IN  THE  SOLID,  WE  HAVE  THE  ADDED  PROBLEM  OF 
AN  UNCOMPENSATED  COUPLE  (ANGULAR  MOMENTUM)  DUE  TO  THE 
DISPLACEMENT  OF  THE  RERADIATED  WAVE. 
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FIGURE  17.  THE  CRITICAL  ANGLE  MEASURING  SYSTEM  CURRENTLY  BEING 

USED  EMPLOYS  A  POINT  DETECTOR  AT  Rj  AND  A  LENS  SOURCE 
AT  S.  IN  THE  MODIFIED  ARRANGEMENT,  WE  INTEND  TO  PLACE 
A  FOCUSED  DETECTOR  (SAME  AS  S)  AT  R2 .  WE  WILL  USE  THIS 
TO  EXAMINE  ANY  ISOTROPIC  RADIATION  COMPONENT  FROM  THE 
REGION  OF  THE  "POINT"  FOCUS. 
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SUMMARY  AND  CONCLUSIONS 


The  Phase  I  study  of  the  "Rayleigh"  critical  angle  phenomenon 
has  been  completed.  An  automated,  computer  controlled  goniometer 
was  constructed  and  perfected.  Controlled  measurements  on  a  variety 
of  materials,  including  materials  of  interest  to  the  NDT  community, 
were  carried  out  and  compared  with  theoretical  predictions  made  by 
a  number  of  authors. 

Experimental  results  indicate  that  these  theories  are  only 
first  approximations  to  reality.  We  find  a  number  of  significant 
experimental  anomalies,  such  as  the  appearance  of  more  than  one 
amplitude  dip.  Certain  theoretical  questions  can  also  be  raised 
regarding  these  theories  including  the  question  of  momentum 
conservation  at  the  interface.  Both  the  experimental  and  theoretical 
observations  suggest  that  we  are  not  yet  in  possession  of  a  complete 
theory . 

Phase  II  plans  are  designed  to  investigate  both  experimentally 
and  theoretically  these  questions  so  that  our  goal  of  a  better 
theoretical  untierstanding  of  this  phenomenon  can  be  achieved.  Given 
this  improved  understanding,  practical  quantitative  NDT  techniques 
based  on  the  "Rayleigh"  critical  angle  phenomenon  will  become 
available. 
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APPENDIX  A:  FUNDAMENTAL  BASIS  FOR  NONLINEAR  NATURE  OF  ELASTIC  WAVES 


Sound  is  an  inherently  nonlinear  phenomenon.  Yet  because  of 

the  possibility  of  expressing  the  motion  of  a  solid  body  in  "normal 

coordinates,"  it  may  be  made  to  appear  as  an  approximately  linear 

phenomenon;  thereby,  obscuring  its  true  nonlinear  origins. 

A  solid  body  may  be  viewed  in  quantum  mechanical  terms 

<*  l 

as  a  collection  of  coupled  harmonic  oscillators.  Each  oscillator 
corresponds  to  an  atom  or  molecule  of  the  solid.  It  is  this  coupling 
that  allows  a  disturbance  at  one  point  in  a  medium  to  propagate  to 
other  points  in  the  medium.  This  coupling  is  essential  for  without 
it  no  energy  or  momentum  could  be  transfered,  and  no  sound  disturbance 
could  propagate.  Such  interactions,  of  course,  render  the  theory 
nonlinear  at  the  outset. 

A  set  of  coupled  quantum  mechanical  oscillators  is  described 

«  1 

by  the  Hamiltonian 
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where  are  the  coordinates  of  the  amount  of  displacement  from 
equilibrium  of  each  oscillator  and  the  momentum  operator  is 
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C.:  are  nonzero  coupling  constants  expressing  the  interactions  among 
the  oscillators.  Making  the  change  of  variables, 
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where  A^  are  constants  independent  of  time,  one  may  ”1 
these  equations  in  the  first  approximation.  One  finds  H 
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where 
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The  eigenvalues  of  the  "collective"  Hamiltonian  H  may  then  be  shown 
to  be 


E 
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(N  +  l/2)-Kw  . 
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This  equation  gives  the  total  energy  of  all  "independent"  a  sound  modes 
present  in  the  solid.  A  particular  sound  mode,  say  a  =  0,  has  energy 


E  =  (N  +  l/2)-hw  . 
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That  is,  the  sound  wave  of  frequency  u)q  consists  of  N0  phonons  each 

with  energy  e  =  -hw  and  momentum 
&  o  o 
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where  v  is  the  phase  velocity  of  the  sound  wave.  In  the  limit  of 


large  quantum  numbers  (N  >  >  1),  the  energy  of  the  sound  wave 
goes  over  to  the  classical  physics  value,  namely 


E  =  N  -hid  =  N  E 
o  o  oo 


and  its  corresponding  momentum  is 


(9) 


(10) 
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which  is  the  usual  classical  expression  for  momentum. 

Strictly  speaking,  if  the  coupling  constants  =  0,  there 
are  no  sound  waves.  Consequently,  they  carry  neither  energy  or 
momentum.  We  see  from  this  exercise  that  although  a  fundamentally 
nonlinear  phenomenon  (sound)  may  look  linear  in  the  classical 
approximation  (in  the  normal  coordinates),  it  is  not.  Thus,  we 
must  use  care,  in  classical  phenomenological  treatments  of  a  sound 
wave,  not  to  eliminate  important  physical  effects  by  using  only 
linear  equations  which  are  fundamentally  incorrect.  The  existence  of 
interactions  between  atoms  in  a  solid,  of  course,  always  leads  to 
inevitable  losses. 

It  is  also  evident  from  the  foregoing  discussion  that  there 
is,  strictly  speaking,  no  lossless  case  in  acoustics.  One  must, 
therefore,  be  careful  in  treating  real  problems  by  first  examining 
lossless,  linear  cases  and  then  making  small  phenomenological 


